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ABSTRACT 

Using VLT/F0RS2 spectroscopy, we have studied the properties of the central (inner 
320 pc) stellar populations of a sample of 26 nucleated early-type dwarf (dE) galax- 
ies in the Virgo Cluster of magnitude range: -18.59 ^ M^. ^ -15.39 mag. With the 
addition of the data of Michiclscn ct al. (2008), extending the sample to 38 dEs, we 
find that these galaxies do not exhibit the same average stellar population charac- 
teristics for different morphological subclasses. The nucleated galaxies without discs, 
which are, on average, fainter than dEs with discs and distributed in regions of higher 
local density, are older and more metal poor (mean ages 7.5 ±1.89 Gyr and 3.1 ±0.83 
Gyr, mean [Z/H] = -0.54 ±0.14 dex and -0.31 ±0.10 dex, respectively). The a- 
element abundance ratio appears consistent with the solar value for both morpholog- 
ical types. Besides a well-defined relation of metallicity and luminosity, we also find a 
clear anti-correlation between age and luminosity. More specifically, there appears to 
be a bimodality: brighter galaxies {Mr ^ —16.5 mag), including the discy ones, exhibit 
significantly younger ages than fainter dEs {Mr ^ —17.0 mag). The magnitude overlap 
between these two subgroups appears to be resolved when considering in addition the 
weak correlation between age and local density, such that older galaxies at a given 
magnitude are located at higher densities. Thus there seems to be no significant differ- 
ence between the stellar populations of dEs with and without discs when compared at 
the same magnitude and density. We thereby revisit the question whether both could 
belong to the same intrinsic population, with discs surviving only in lower-density 
regions. However, it appears less likely that fainter and brighter dEs have experienced 
the same evolutionary history, as the well-established trend of decreasing average stel- 
lar age when going from the most luminous ellipticals towards low-luminosity Es and 
bright dEs is broken here. The older and more metal-poor dEs could have had an 
early termination of star formation activity, possibly being "primordial" galaxies in 
the sense that they have formed along with the protocluster or experienced very early 
infall. By contrast, the younger and relatively metal-rich brighter dEs, most of which 
have discs, might have undergone structural transformation of infalling disc galaxies. 

Key words: galaxies: dwarf - galaxies: evolution - galaxies: formation - galaxies: 
stellar content - galaxies: elliptical and lenticular, cD - galaxies: clusters: individual: 
Virgo 



1 INTRODUCTION 

Dwarf elliptical (dE) galaxies are the nu merically domi- 
nant population in the present-day Universe jSandage et al.l 



* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile 
(programme 078.B-0178) 
f E-mail: sjy@x-astro.net 



Il985l : iBinggeli et allTigSTh . They also exhibit strong clus- 
tering, being found predominantly in the close vicinity 
of giant galaxies, either as satellites of individual giants, 
or a s members of galaxy clust ers (|Ferguson fc Sandagd 
1 19891 : iFerguson fc Binggelil Il994l '). Unlike classical ellipti- 
cal galaxies, stellar population studies show that these 
galaxies exhibit on average younger ages as compared to 
their giant counterparts, and also a lower metal content 
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acco rdinR to the correlati on of metallicity and luminos- 
ity (|Michielsen et all |2008| . hereafter M08). More detailed 
studies of dEs and the fainter dSphs in the Local Group, 
based on resolved stellar populations, reveal that most 
dE/d Sphs have a fairly extended star for mation history 
(e.g., lMateolll998l : iGrebel fc GallagheJiiooi ). with the last 
star formation activity ranging from a few gigayears ago 
to the time of re-i onisation. A number of studies in diffe r- 
ent clusters (e.g.. iPoggianti et al]|200ll: iRakos et al.l[200 



stripping jGunn & GottI 


19721: Faber & Lin 19831: Silk et all 


ll987l:lBabul & Reeslll992 


:lMurakami & Babulll 19991) and are 



ICaldwell et~alll2003l : iGeha et al.ll2003l : Ivan Zee et al.([20oi ) 

also provide a wide range of ages for passive dwarfs. An- 
other study, h owever, based on the velocity distribution of 
the Virgo dEs (|ConseIice et al ] |200lh . argue that dEs are not 
an old cluster population. This evidence suggests that dEs 
in nearby clusters have a mixed origin. Some have proper- 
ties consistent with an early "primordial" formation, while 
others appear to be more recently formed, possibly from the 
transformation of progenitors with different morphological 
types. 

Due to the distance of the Virgo cluster (16.5 Mpc, 
iMei et cd] |2007| ). it is basically impossible to perform re- 
solved stellar population studies of dEs in Vi rgo (see, 
howe ver, the studies of the red giant branch by ICaldwelll 
l2006h . Therefore, in order to probe the star formation his- 
tory of these faint galaxies, we use integrated spectra with 
a good signal-to- noise ratio (SNR), providing a sensitive 
tool for stellar population studies through absorption fea- 
ture measure ments. The basic id ea, which was first in- 
troduced by iBurstein et al.l Il984l . is to provide a set of 
optical absorption line strength indices, co mmonly known 
as the Lick/IDS system or Lick indices ([Burstein et al.l 
1 1984 IWorthe7 et al. l ll994al : [TVager" et al.l 0*9981 ). It is now a 
widely used approach to determine the age, metallicity, and 
Q— element abundance of galaxies from integrated spectra. 
The measured indices from the optical spectrum of a galaxy 
are compared to their values predicted by stella r popula- 



tion models, which are provided for example by IWortheyl 
1994 IVazdekis fc Arimotd Il999l : iBruzual fc CharlotI l200a 



Thomas et al .ll2003l : ISchiavonl2007l ~ The SSP-equivalent age 
and metallicity of a given galaxy are then determined as the 
ones of the model whose indices best agree with the mea- 
sured ones. 

It has already been noticed that dwarf elliptical galax- 
ies of the Virgo c luster are consiste nt with solar [a/Fe]- 
abundance ratios (|Gorgas et al.lll997l ). indicating the slow 
chemical evolution in low-mass sy stems. Other studies 
l|Geha et al.ll2003l : lvan Zee et al.ll2004l ) confirmed this result. 

As the most numerous type of galaxy in clusters, 
early-type dwarf galaxies are ideal probes to study the 
physical processes that govern galaxy formation and evo- 
lution in environments of different dens ity. The pro- 
nounced morpholog y-density relation (e.g.. IPressleil [l980l : 
iBinggeli et al.l 119871 ) suggests that early-type dwarfs were 
either formed mainly in high-density environments, or orig- 
inate from galaxies that fell into a cluster and were mor- 
phologically transformed. Because having low masses their 
properties are expected to depend strongly on the environ- 
ment they reside in. However, the actua l formation mech- 
anism s are still a matter of debate (see Ijerien fc Binggelil 
I2OO5I . and references therein). Most of the proposed scenar- 
ios are based on the vigorous forces acting within a cluster 
environment, such as interactions with a hot intracluster 



able to transform dw arf irregular (dirr) g alaxies or late-type 
spirals into dEs fe.g.. lvan Zee et al ]|2004i. Another scenario, 
the so-called harassment ( Moore et al.lll99^ . can transform 



infalling late- type spirals through close e ncounters wi t h mas - 
sive cluster members , and the studies of Jerien et al.l (|200Ci ) , 
iBarazza et al] (|2002l ). and lDe Riicke et al.1 (|2003l ) also sug- 



gest the possibility of production of dEs by harassment. 

A recent systematic study bv lLisker et al.l (|2006al . [20071 . 
hereafter LOB and L07, respectively) revealed a striking het- 
erogeneity of the class of early-type dwarfs. They found dif- 
ferent subclasses, with significantly different shapes, colours, 
and spatial distributions. Those dEs with a disc component 
have a flat shape and are predominantly found in the out- 
skirts of the cluster, suggesting that they - or their progeni- 
tors - might have just recently fallen into the cluster environ- 
ment. In contrast, dEs with a compact stellar nucleus follow 
the classical picture of dwarf ellipticals: they are spheroidal 
objects that are preferentially found in the dense cluster cen- 
ter. Thus, given the structural heterogeneity and the scatter 
in their stellar population characteristics, the question arises 
whether they have a similar or different formation history? 
Here we address this question with exploring the relation be- 
tween their stellar populations, structural parameters, and 
the environment. 

The paper is structured as follows. In Section 2, we 
provide an overview of the data, the sample selection, and 
the observational parameters. In Section 3, we present the 
basic data reduction process. We describe in Section 4 the 
steps followed towards the Lick index measurement, includ- 
ing consistency checks of the models used, as well as tests 
and comparisons of Lick indices at different spectral reso- 
lutions. The derivation of stellar population parameters is 
outlined in Section 5. The results from our analysis are pre- 
sented in Section 6, separating the actual index measure- 
ments (Section 6.1) and the subsequently derived ages and 
metallicties (Section 6.2). Finally, Section 7 presents the dis- 
cussion of the results, and our conclusions. 



2 SAMPLE AND OBSERVATIONS 
2.1 Sample selection 

Among the cert ain cluster members in the Virgo cluster 
catalogue fVCC. [Binggeli et al.|[l985l : IBinggeli fc CameronI 
Il993h . we selected nucleatecf] dEs (dE(N)s) of different 
morphology, according to L06. Our target list comprises 8 
dEs with discs (dE(di)s), 7 dE(N)s without discs located 
in the central region of the cluster, and 11 dE(N)s with- 
out discs located in the outer cluster regions (including 4 
objects in the dense southern subcluster). The majority of 



^ The classification of a dE as nucleated and non-nucleated is 
not unambiguous. Many Virgo dEs classified as non-nucleated in 
the VCC ac tually host a faint nucleus hardly de tectable with the 
VCC data llGrant et al.ll2005l : ICote et al.ll2006ll A more appro- 
priate term might be dEs without a nucleus of significant relative 
brightness as compared to the central galaxy light. 
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Figure 1. Colour-magnitude diagram of Virgo cluster dEs with 
discs (green), without discs (red) and blue centre (blue) as iden- 
tified by L07. Solid circle represent the dEs of our target sample. 



dE(di)s populate the outer cluster regions (L07) and there- 
fore need to be compared to similarly located dE(N)s not 
hosting a disc. One of the dE(di)s of our sample is known to 
host young stars in its centra l region, reflected in the blue 
central colour (VCC 0308, see lLisker et al.ll2006bh . 

One major goal of our observations was to be able to 
analyse the nuclei themselves, which will be subject of a 
second publication to follow. Therefore, our target selection 
was guided by the estimated ratio of nucleus light to un- 
derlying galaxy light, and by the SNR of nucleus and host 
galaxy. Our sample is therefore mainly biased towards galax- 
ies with relatively bright nuclei, both in absolute terms and 
relative to the galaxy's central region. We will take this into 
account in the discussion section. 

Furthermore, as shown in L07, the dE(di)s have, on 
average, somewhat bluer colours than the dE(N)s. This is 
reflected in our sample: Fig. [T] shows that the dE(di)s tend 
towards the bluer side of the colour-magnitude relation in 
g-i. Whether this effect is due to age or metallicity, and 
whether it is caused by the presence of a disc itself or by 
environmental density, will be addressed in the discussion 
section. 

We believe that this is a representative sample of bright 
nucleated early-type dwarf galaxies of the Virgo cluster with 
different substructural types. Table[T]lists the sample galax- 
ies and some of their basic properties. We also build an 
extended sample by adding Virgo dEs from the M08 sam- 
ple whose morphology is defined in L06, and remeasuring 
their stellar population parameters from the reduced spec- 
tra, which were kindly provided by C. Conselice on behalf of 
the MAGPOP ITF0 collaboration. M08 selected their galax- 



ies to be brighter than ms <~ 15 mag, with high central 
surface brightness, and all but one (VCC 0917) hosting a 
central nucleus. In total, our extended sample thus com- 
prises three dE(bc)s (of which one also hosts a disc, as men- 
tioned above), 13 dE(di)s, and 22 dE(N)s that do not show 
disc features. Note that we removed VCC 0917 from M08 
sample because it has a different subclass, as this galax;y is 
neither discy, nucleated, nor blue-centered according to L06. 



2.2 Observational characteristics 

The observations were carried out over six half nights 
from March 16 to 22, 20 07, with ESO VLT UT1/F0RS2 
(|Appenzeller et al] 1 19981 ) in multi object spectroscopy 
(MXU) mode. The dates and other observational informa- 
tion are provided in Table [21 In this table, the last two 
columns give the number of exposures, the exposure time, 
and the resulting SNR per pixel of the combined spec- 
tra, measured at 5000 A rest frame wavelength after co- 
adding them. Spectrophotometric standard stars were ob- 
served each night. The typical value of the seeing FWHM 
during the galaxy exposures was 1.3" as measured by Gaus- 
sian fitting of the field stars. The detector in F0RS2 com- 
prises two 2k X 4k MIT CCDs, with a pixel size of 15 
X 15 fj,m^. The standard resolution collimator and the 2- 
pixel binned read-out yield a image scale of 0.25"pixel^^ . 
The spectrograph slit width was chosen to be 1", and cov- 
ered 40" in length, using the GRIS300V grism providing 
a dispersion of 3.36 Apixel"^. This setup yields a spec- 
tral resolution, as measured from the FWHM of the arc 
lines, of ~11 A at ~5000 A, which is below the instrumen- 
tal resolution of Lick/IDS system (~8.4 A at 5000 A) see 
IWorthev fc Ottaviani|[l997l . 



3 DATA REDUCTION 

The first step in the data reduction was performed with 
the VLT F0RS2 Pipeline GasgancO implemented on the 
EsoRe>0 data reduction package. This software performs 
the following basic tasks: bias subtraction, aperture identifi- 
cation, flat-fielding, cosmic-ray removal, and wavelength cal- 
ibration (using the arc lines of the HeNe lamp). The recipe 
fors_caUb is used for the definition of the extraction mask 
based on the flat field and arc lamp exposures, and is also 
used for the creation of the normalized flat and master bias. 
In addition to that, fors.caUb also calculates the dispersion 
coefficients which are used for wavelength calibration during 
the final reduction of the science frames using the second 
recipe fors.science. The typical uncertainties of the wave- 
length solution were around 0.2 pixels or less, which is a 
good agreement with the FORS pipeline standard. 

Since our spectra were taken in MXU mode, which al- 
lows a maximum slit length of 40" , the slit does not cover the 
full extent of our target galaxies and a sufficient region of the 
sky at the same time; instead, additional slits were placed 
on blank sky regions. For this reason, the F0RS2 Pipeline 



2 http://www.astro.rug.nl/~peletier/MAGPOP.ITP.html 



^ http:/ /www. eso.org/sci/data-proeessing/software/gasgano 

* http://www.eso.org/sci/data-processing/softwarc/cpl/esorex.html 
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Table 1. The sample: basic properties. 



Name RA Dec Type Density_^ K^g RV Remark 

1 2 3456~7 8 9 
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21.0 


22.5 


1871 ± 16 


also in M08 


VCC1304 


12: 


:30:39.90 


15:07:46.68 


d,n 


-16.86 


15.1 


16.2 


—109 it 78* 




VCC1308 


12: 


:30:45.94 


11:20:35.52 


n 


-16.50 
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extended 
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VCC0021 


12: 


:10:23.15 


10:11:19.04 


b,n 


-17.01 


1.8 


15.2 


486 ± 25 




VCC0397 


12: 


:20:12.18 


06:37:23.56 


d,n 


-16.76 




12.4 


2471 ± 46 


poss. mem. 


VCC0523 


12: 


:22:04.14 


12:47:14.60 


d,n 


-18.57 


20.2 


25.7 


1981 ± 20 




VCC0917 


12: 


:26:32.39 


13:34:43.54 




-16.55 


45.9 


9.2 


1238 ± 20 


removed 


VCC1087 


12: 


:28:14.90 


11:47:23.58 


n 


-18.50 


54.9 


28.4 


675 ± 12 




VCC1122 


12: 


:28:41.71 


12:54:57.08 


n 


-17.13 


66.0 


19.1 


476 ± 10 




VCC1183 


12: 


:29:22.51 


11:26:01.73 


d,n 


-17.82 


28.2 


19.6 


1335 ± 12 




VCC1431 


12: 


:32:23.41 


11:15:46.94 


n 


-17.72 


30.2 


9.3 


1505 ± 21 




VCC1549 


12: 
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11:04:17.51 


n 


-17.24 


28.3 


10.9 


1357 ± 37 
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21.4 
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1547 ± 29 
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d,n 


-17.86 


28.1 


13.6 


206 ± 26 
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12: 
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b,n 


-17.83 


19.6 


22.3 


-169 ± 28 




VCC1947 


12: 


:42:56.37 


03:40:36.12 


d,n 


-17.61 




9.4 


974 ± 19 


poss. mem. 



In column 1 we give the galaxy name according to the VCC. Column 2 and 3 represent the galaxy position 
in RA and DEC (J2000.0), respectively. For column 4 we take the morphological classification given by L07 (with 
b=blue-centre, d=discy n=nucle ated). Column 5 gives the galaxy absolute magnitude in the SDSS r-band taking 
m— M = 31.09 l lMei et al.ll200'^ . corrected for Galactic extinction (L07). Column 6 and 7 represent the local 
projected density (number per sq.deg., calculated from a circular projected area enclosing the 10th neighbor), 
and the efi^ective radius in arcsec measured by L07, respectively. In column 8 we give the radial velocity in km/s 
provided by NED; for those galaxies whose radial velocity is not listed in NED* we measured their radial velocity 
from our data set, albeit with larger errors due to the low spectral resolution. In the last column we provide 
comments and indicate galaxies in common with the study of M08. 



Local projected galaxy number density 



does not provide a good result for the sky background sub- 
traction for this type of spectra of spatially extended objects 
(as also mentioned in the pipeline manual), and left signifi- 
cant residual signatures of sky lines. Therefore, the night sky 
line removal has been done with iRAljfl using the background 

^ Image Reduction & Analysis Facility Software distributed by 



task implemented on the oned package. Inherently, residuals 
and uncertainties from the sky subtraction could affect the 

National Optical Astronomy Observatories, which arc operated 
by the Association of Universities for Research in Astronomy, Inc., 
under co-operative agreement with the National Science Founda- 
tion 
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Table 2. Observation Log 



Name 


Chip 


Night 


Exp. (s) 


SNR(pixel-i) 


VCC0216 


1 


2007-03-19 


2 


X 


660 


46 


VCC0308 


1 


2007-03-22 


6 


X 


590 


49 


VCC0389 


1 


2007-03-22 


2 


X 


510 


48 


VCC0490 


1 


2007-03-18 


3 


X 


420 


36 


VCC0545 


1 


2007-03-17 


3 


X 


680 


43 


VCC0725 


1 


2007-03-18 


6 


X 


670 


31 


VCC0856 


2 


2007-03-19 


1 


X 


780 


58 


VCC0929 


1 


2007-03-19 


1 


X 


640 


63 


VCC0990 


1 


2007-03-22 


3 


X 


660 


61 


VCC1167 


1 


2007-03-18 


3 


X 


720 


45 


VCC1185 


1 


2007-03-20 


2 


X 


705 


42 
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2007-03-18 


2 


X 


600 


56 


VCC1261 


2 


2007-03-20 


1 


X 


840 
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2007-03-17 
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X 


570 


43 
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2007-03-20 
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X 


540 


41 
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2007-03-17 


4 


X 


570 


38 
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2007-03-19 


2 


X 


450 


43 
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1 


2007-03-19 


6 


X 


760 


40 
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1 


2007-03-21 


5 


X 


670 


34 
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1 


2007-03-20 


6 


X 


670 


40 
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2 


2007-03-22 


5 


X 


680 


46 


VCC1661 


1 
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dard star spectrum and the tabulated calibrated spectrum, 
and the resulting continuum correction was subsequently ap- 
plied to the galaxy spectra. To correct for the effect of at- 
mospheric extinction, we used the iraf task setairmass to 
define the effective airmass for each spectrum. Then, we ap- 
plied the extinction correction (task calibrate) using the ap- 
propriate atmospheric extinction table (provided by Sabine 
Moehler) . The spectral sensitivity functions as derived from 
the different stars varied by typically reaching ~3% to- 

wards the edges of the wavelength range. Given these small 
variations, we created a single function from all standard 
stars, and are reasonably confident that the resulting inter- 
nal relative flux calibration is good to 

The individual spectra from different exposures of a 
given galaxy were then co-added using the iraf task scom- 
bine to produce a spectrum of higher SNR. Three typical 
flux calibrated combined galaxy spectra are shown in Fig. 

m 

For the later calculation of index uncertainties, we 
also produced error frames, following the description by 
ICardiel et al.l (|l998l ). in particular the processing of error 
frames parallel to the reduction steps. In this sense, we have 
an error spectrum for each fully processed science spec- 
trum that contains the propagation of the initial random 
errors (due to photon statistics and read-out noise) carried 
throughout the arithmetic manipulations in the reduction 
procedure. 



line strength measurements that we will obtain from these 
data, so we have investigated thoroughly how this step is im- 
plemented. The sky background is obtained from the blank 
slits, which were placed along the extension of the galaxy 
slit. All the exposures have at least one blank slit for the 
sky background, in most cases more than one. Nevertheless, 
we found that the resulting sky-subtracted galaxy spectrum 
was almost completely insensitive to how the sky region was 
chosen, and whether one or two slits were used to sample 
the background. We have therefore chosen the sky regions 
using the full length of the blank slit, typically 38" , omitting 
1" at either side of slit. 

The photometric standard stars were observed in MOS 
mode. Here, the FORS2 pipeline sky subtraction did not 
leave any residuals, since the target is a point source and 
a sufficient amount of blank pixels are thus covered by the 
slit. 

Once the galaxy frames were fully reduced, one- 
dimensional spectra were extracted by summing the pixel 
values in spatial direction along the central four arcsec- 
onds, corresponding to 3 20 pc with a d istance modulus 
m- M = 31.09 (16.5 Mpc. lMei eralll2007l 'l . This is done in 
order to guarantee that the resulting Lick indices will have 
sufficiently small errors to allow a precise stellar population 
analysis, and to enable a direct comparison with the sample 
of MOS, who used the same extraction aperture. This allows 
us to build an extended sample (see Sect. I2.l"|) by combining 
their and our galaxies. 

The spectrophotometric calibration was done using 
the observations of 5 standard stars (LTT4816, LTT6248, 
LTT1215, LTT2415, and Feige56). The instrument response 
function was computed from the ratio of the observed stan- 



4 LINE INDICES: MEASUREMENT AND 
CALIBRATION 

In this section we provide the details of the measurement 
of line strength indices in the wavelength range 4000 A to 
6000 A. The main step is to calc ulate the line stre ngths 
particularly for Lick/IDS indices (|Trager et al.l [l998l ). us- 
ing one-dimensional spectra. This will allow us to com- 
pare the results t o simp le stellar populat ion models of 
iBruzual fc CharlotI l|2003l hereafter BC03), iThomas et al] 
(|2003l . here after TMB03) and AV/MILES (Vazdekis et al. 
2009; in preparationjf], in order to derive SSP-equivalent 
ages and metallicities. In this work, we focus on maximiz- 
ing the SNR of the spectra and on the ability to combine 
our dataset of 26 galaxies with the 12 of M08, who obtained 
their 1-D spectra by summing the central 4 arcsec in spa- 
tial direction. We therefore apply the same extraction to 
our spectra, and defer any stu dy of potenti al population 
gradients to a later publication. iLisker et al.l stated that a 
definite conclusion, as to whether the observed colour differ- 
ences of different morphological subtypes are caused by age 
or metallicity effects, was not possible. With our good SNR 
in the central region, we can address this issue in a more 
appropriate and reliable way, by comparing our index mea- 
surements with stellar population models, in order to break 
the degeneracy of age and metallicity. 
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Figure 2. Flux calibrated spectra of three typical galaxies in our sample are shown. The specatra are shifted to restframe wavelength 
and normalized at 5500 A. Selected, important Lick/IDS indices central band pass are indicated by the gray shaded region. 



4.1 Line strength Measurement 

The hne index strengths were measured from flux cahbrated 
spectra using the routine Indexi{3 developed b y N. Cardiel. 
It us es the definition of the Lick indices from Trager et al.l 
Il998l and calculates a pseudo-continuum (a local continuum 
level) for each spectral feature defined by the mean values 
within two pseudocontinuum bandpasses on either side of 
the spectral feature. This software also estimates the un- 
certainties resulting from the propagation of random errors 
and from the effect of uncertainties on radial velocity by 
performing Monte-Carlo simulations. 

Note that we have not applied a velocity dispersion cor- 
rection for our galaxies, because the expected galactic veloc- 
ity dispersion, agai ^ 50 kms~^, is significantly below our 
spectral resolution ainstr ~ 280 kms^^. Therefore these cor- 
rections are not necessary. 



4.2 Index Calibration and Stellar Population 
Model 

There are several stellar population models that provide Lick 
indices fo r stellar popula t ions of different ages and m etallic- 
ities (e.g. IWorthe"^ll994l : IVazdekis fc Arimotalll999l : BC03; 
TMB03). The models adopted in this paper are those of 
TMB03, BC03, and Vazdekis/MILES (new version based 
on MILES library, hereafter AV/MILES). The model of 
TMB03 has a variable a-element abundance ratio, based 
on the st a ndard SSP m odels computed with the code of 
iMarastonI (| 19981 ) with a ISalpeteij (|l955h IMF. They pro- 
vide Lick index values at Lick standard resolution. The 
new AV/MILES models are ba sed o n the previous model s 



19991) andlVazdekis et all(|2003h . 
Girardi et all I2OO0I ') isochrones. 



from lVazdekis fc Arimotol 
They use Padova 2000 
They cover ages between 0.1 and 17.8 Gyr, and metal- 
licities between -1.7 and .2 dex. The MILES library 
(|Sanchez-Blazauez et al.ll2006l ) used for these models has a 
2.3 A spectral resolution (FWHM). It covers the wavelength 
range from 3525 to 7500 A, with a good calibration of the 
flux. 

Among these three models (i.e., BC03, TMB03 & 



http:/ /www. ucm.es/info/Astrof/software/indexf/indexf.html 



AV/MILES), we used BC03 and AV/MILES for the con- 
sistency test (see next paragraph), since they provide the 
full spectral energy distribution (SED) of SSPs, and can 
therefore be easily transformed to our spectral resolution 
of 11 A. Furthermore, we also use AV/MILES for the cor- 
rection of measured indices back to the Lick system (see 
section iLSl for details of the correction procedure). We use 
the x^— technique to derive stellar population estimates from 
AV/MILES at 11 A resolution and TMB03 at Lick resolu- 
tion (see section [5]). Finally, after conflrming the consistency 
between the stellar parameter estimates from these two dif- 
ferent models, we use those obtained with TMB03 as our 
flnal values, since the TMB03 model also covers different 
[a/Fe] ratios, and thus avoids potential problems of other 
models with flxed solar abundance ratio. 

Before showing the main results, we provide consistency 
tests between our data and the different SSP model predic- 
tions of BC03 and AV/MILES. For this purpose, the model 
indices were computed at the resolution of our data (~11 A). 
We then create index-index plots that c ontain those indice s 
sensitive to the same chemical species (|Kuntschneill2000l '). 
These plots provide an idea of how accurately the models 
reflect the data and, in consequence, they help us to inter- 
pret the results with respect to a chosen model. Moreover, it 
is important to check whether or not signiflcant differences 
exist in the relative spectrophotometric calibrations of the 
data and the stellar libraries used by the models. 

In Figure [3] different Mg indices are plotted against 
each other. The folded grids represent BC03 (black) and 
AV/MILES (red) models at 11 A resolution, and the black 
dots are the indices measured from our sample. The plots 
show a good overall agreement between the measured index 
values and the model predictions, with a small systematic 
deviation for galaxies with smaller Mg indices. Likewise, Fig. 
Ushows the relation between <Fe>, Fe5015 and Fe5406. In 
this plot, the scatter of the data points around the model 
grid locus is larger; however, this is explained by our mea- 
surement errors being larger there. Overall, the models do 
trace the mean trends in the data as a whole. In Fig. [5] we 
present the Balmer indices II/3, II5_f and II7F. Here we flnd 
fairly good agreement for the AV/MILES model, while a 
small deviation can be seen between the BC03 model pre- 
dictions and the data in the HiJf and II7F indices. 
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4.3 Correction to the Lick system 

In this section, we investigate tlie beliaviour of line indices 
at different resolutions, and provide corrections of our mea- 
sured indices to the Lick system. We assess possible effects 
of a variation of the resolution on the ages and metallicities. 
For this purpose, we used AV/MILES model spectra. First, 
we degraded the model spectra to Lick resolution, since most 
studies in the literature are based on this resolution. How- 
ever, for simplicity, we consider a fixed resolution of 9 A, 
although the original Lick system varies from ~8.4 to 11.5 
A depending on wavelength. In addition, we computed an- 
other degraded version of the model spectra, adopting the 
resolution of our data (11 A). These two different resolu- 
tions are compared with each other in Fig|S] where we show 
the difference AI between the measured indices for the two 
resolutions (i.e., AI = Indices measured at 9 A — 11 A) as a 
function of index strength. We find that Fe5335 and Fe5406 
suffer slightly more degradation in the strong absorption re- 
gion than the Balmer lines (H/3, H7F, and H5_f) and the 
Mg indices. The percentage variations of these features are 
-10.5% for Fe5335 & Fe5406, ~6% in case of Fe5270 & 
Fe5015 and Mgb, and even smaller for II/3 (i.e., ~2%) and 
other higher order Balmer lines. Note that these system- 
atic deviations in absorption line strengths are smaller than 
half of the typical (statistical) measurement error from our 
galaxy spectra. 

From the relations illustrated in Fig|B] we computed 
corrections for index strength for all the measured indices 
of our galaxy sample, in order to correct our data back to 



the Lick resolution (9 A) . This allows us to apply a separate 
correction for each galaxy, depending on its individual index 
strength. We then further corrected each index by applying 
the offsets between flux-calibrated spectra and the non-flux- 
calibrated original Lick/IDS system. For this purpose, we 
derived offsets from the Lick star spectra of the MILES li- 
brary (see the Section|X]for more detail). With these correc- 
tions, our data qualify for comparisons with various models 
based on the Lick resolution. We are thus in the position 
to compare our data to the TMB03 models, which provide 
tabulated values of Lick indices as a function of age, [Z/H] 
and abundance ratio [a/Fe]. 

To examine the robustness of our approach of correcting 
the indices back to the Lick system, we comp are the final 
corrected values of those galaxies common to iGeha et al.l 
(200J) and M08 with their published values. In Fig. [3 we 
find that these adopted corrections give results consistent 
with the previously published data. Therefore, our further 
analyses (i.e., index analysis, comparison to models, and ex- 
traction of age and metallicity) are done on the Lick system, 
unless explicitely stated otherwise. 



5 EXTRACTION OF STELLAR POPULATION 
PARAMETERS FROM MODEL FITTING 

As the aim of our work is to understand the star formation 
and evolutionary history of dwarf galaxies, we extract stellar 
population characteristics from those indices, and combina- 
tions of indices, that are sensitive to age, metallicity, or the 
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Figure 4. Consistency test using Fe indices and BC03 models at 11 A resolution. The folded grids are BC03 (black) and AV/MILES 
(red) models and the datapoints correspond to the central indices measured in our sample. 




Index strength at 11 A 



Figure 6. The deviation of index values after changing the reso- 
lution from qA to llA as a function of measured index strength 
at llA of AV/MILES model spectra. 

relative abundance of different metals. Population synthe- 
sis models then provide us with SSP-equivalent, or mean 
luminosity- weighted, ages and metallicities. It is necessary 
to mention th at the SSP-equivalent ag es are biased towards 
younger ages (|Trager et al.ll2000l , lioosl) since in the optical, 



most of the hght comes from the youngest component of 
the stellar population, and the age correlates most strongly 
with last star formation activity. Therefore, the measured 
age is influenced to a large extent by the last star formation 
activity in the galaxy. However, the SSP-equivalent metal- 
licity is an excellent tracer of the light-weighted metallicity 
PTr ager fc Somervillell2009l ). This is because hot, young stars 
contribute little to the metal lines in a co mposite spectrum 
(jTrager et al.ll2005l : ISerra fc TrageilbOOTh . 

To convert our measured indices to ages and metallici- 
ties, we used the y ^— minimization method as suggested by 
iProctor et al ] |2004l . In this method, we have the advantage 
to use a large set of indices to get the best fitting stellar pop- 
ulation parameters. Therefore, it makes better use of the 
information that directly depends on age, metallicity and 
a— element abundance ratio, as compared to an analysis of 
single indices or index pairs. To perform the model fits, we 
first interpolated the model grid of TMB03 to a finer grid. 
The actual x'^— minimization was then performed on the age, 
[Z/H], and [o/Fe] parameter space. While the method al- 
lows us to use as many indices as possible, we used only 
the nine indices with the best measurement quality (H5f, 
H7F, Fe4383, H/3, Fe5015, Mgb, Fe5270, Fe5335 & Fe5406) 
among the measured Lick indices from our spectra, which 
are tabulated in Table [Bl] To calculate the uncertainties in 
the stellar population parameters, we transformed the index 
error to the stellar population error contours and obtained 
the maximum of the error from that. 
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Figure 5. Consistency test using Balmer indices indices and BC03 models at 11 A resolution. The folded grids are BC03 (black) and 
AV/MILES (red) models and the datapoints correspond to the central indices measured in our sample. 




In Fig. [8] we show the comparison of derived ages 
and metallicities using different SSP models (TMB03 & 
AV/MILES). We derived TMB03 ages and metallicities by 
correcting our data to the Lick system as described above, 
while AV/MILES ages and metallicities were obtained di- 
rectly at the resolution of our data (11 A), by degrading 
the AV/MILES model down to this resolution prior to the 
model fitting. We use the same set of nine indices in both 
cases. Although, the measured age and metallicity agree well 



within the error limit, it seems that the majority of derived 
ages are overestimated, when we used the model TMBo|f|. 

* As a further test to check the source of this slight offset, we com- 
pared the estimated SSP parameters from the model AV/MILES 
at different resolutions (i.e., at 9A and llA) for the M08 sam- 
ple. While no strong deviation is seen, there might be a small 
systematic bias of having younger ages when we degrade the res- 
olution (the plot is given in Appendix [C] Fig.[C4ll. In any case, 
the estimated ages at different resolution agree well within the 
errors. 
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The behavior of the nine well-measured indices with re- 
spect to the best-fit model index value is shown in Fig. (9] 
Each colour represents a different index; the galaxy number 
along with its subtype label is provided as y-axis. Clearly, 
the index values lie well within the 2a uncertainty range, 
and surprisingly, in the vast majority of cases even within 
the 1(7 limit, ft is interesting that the deviation of the Mgb 
indices for the AV/MILES model is systematically negative 
in almost all cases, and only for VCC 1348 it shows a sig- 
nificantly positive deviation. This galaxy has the highest 
a— element abundance among our sample (see Table I13p . 
On the other hand, for the TMB03 model, the distribution 
of Mgb indices is rather symmetric, and all lie within the la 
error limit. Therefore the asymmetric distribution of Mgb 
for the AV/MILES model can be interpreted as an effect of 
having a fixed solar a— element abundance ratio. We also 
note that the deviations of Fe5015 and Fe4383 are mostly 
positive, while Fe5335 and Fe5406 have negative deviations 
in most cases. However, the overall mean of the deviations 
for all indices is always less than 0.5 sigma for the model 
TMB03. 

The derived values of age and metallicity for our dEs, 
along with their uncertainties, using the TMB03 model, are 
tabulated in Tabled We also remeasured age and metallicity 
for the spectral sample of M08, using the same nine indices 
as in our sample. This guarantees a proper comparison, be- 
cause M08 have used a different approach (i.e. quadratic 
interpolation over the nearest SSP model grid points in 
the H/3-[MgFe] dia gram, as suggested by the method of 
ICardiel et al.|[200^ ). Interestingly, we obtain a smaller error 
on age and metallicity, as compared to their measurement. 
Nevertheless, the measured ages and metallicities match 
very well (see Fig. IClf) . 



6 RESULTS 

In this section, we first show the measured indices directly, 
in correlations with galaxy absolute magnitude, and in sim- 



ple diagnostic index vs index plots, which compare them to 
model grids of age and metallicity. Finally, we analyse the 
stellar population results obtained from model fitting. 

6.1 Lick indices for Virgo dEs 

The relation of the galaxy magnitude and the in- 
dices that are widely used as age and metallicity in- 
dicators (i.e. H/3 as good age indicator and [MgFe]' 
as good mean metallicity indicator, wit h [MgFe]^ = 
x/M^h X (0.72 X Fe5270 + 0.28 x Fe5335) (|Thomas et al.1 
(2OO30) are shown in Figure [10] It is remarkable that we do 
not see any correlation between the age sensitive index H/3 
and magnitude, while a small but clear offset between nucle- 
ated early-type dwarfs with and without discs (dE(di) and 
dE(N), respectively) is seen. On the other hand, the metal- 
licity sensitive index [MgFe]' shows a weak anti-correlation 
with Mr, from the b righter dE(di)s down to the fainter 
dE(N)s. As we know (jWorthev et al.l [l994bl ') . there are no 
such pure indices which only depend on either age or metal- 
licity. Therefore, the exploration of age and metallicity by 
comparing the indices or their combinations to model pre- 
dictions is presented further below. 

VCC 0725 is an outlier, due to its fairly low Mgb content 
(see Table A2). This galaxy is one of the faintest among 
our sample and there is almost no galaxy light beyond the 
central nucleus. Therefore, we suspect that a fairly large 
domination of the residual sky noise within the four arcsec 
aperture could produce such an effect. Hence, we remove 
this galaxy from the sample during the subsequent analysis. 

Figure [TT] shows the distribution in H/3 versus [MgFe]' 
of dEs and a grid of simple stellar population models of 
BC03, AV/MILES & TMB03 (with [o/Fe] = 0.0 dex). We 
use [MgFe]' as a metallicity indicator, which is considered 
independent of [a/Fe]-abundance, and H/3 as age indicator, 
because its age sensitivity is greater and it is less degenerate 
with metallicity or abundance rat io variations than the other 
Balmer indices (|Korn et al.ll2005l ') . 

The solid lines, which are close to vertical, are of con- 
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Figure 9. Distribution of tlie deviation of indices from tlie best fit model point. The Y-axis is labeled with the VCC number of our 
target galaxies preceded by a letter indicating their galaxy type. 
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Figure 11. The ago-sensitive index H/3 as a function of the metallicity-sensitive index [MgFe]'. Overplotted are the stellar population 
models of BC03(left), AV/MILES (middle) & TMB03 with [a/Fc]=0.0 dex (right). The points without error bars are from the M08 
sample; the mean error for these data is shown as a cross in the upper right corner of each panel. Colours are the same as in Fig llOl 



stant metallicity, whereas the dashed lines, which are close 
to horizontal, connect constant-age models (from top to bot- 
tom: 1, 2, 3, 5, 7, 10, 12, & 15 Gyr). All three panels appear 
to be fairly consistent with each other. However, there are 
some galajcies that lie outside of the model grid, located at 
low H/3, i.e., at the bottom left corner of the grid. Neverthe- 
less, when taking into account their errorbars, they could 
still lie within the grid region. Moreover, it is not surprising 



that some galaxies (particularly dEs) fall outside of model 
grid. The same problem has already been noticed in the 
similar type of study of the dEs in the Coma cluster by 
IPoggianti et all (|200ll '). This can be an issue of how the mass 
loss along th e red giant branch is treated in the models, as 
discussed bv lMaraston et aDl2003l . 

We can see that most of the dEs appea r to have sub- 
solar metallicity, supporting the results of Ivan Zee et al.l 
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Figure 10. Relations between Mr and the indices used as age and 
metallicity indicators (dE(bc): blue, dE(di): green, dE(N): red). 
Open circles represent dEs from the M08 sample. The median of 
the indices is represented by the dashed horizontal line for the 
dE(N)s and by the solid line for the dE(di)s. 



l|2004 ). M08, iGeha et all l|2003h . iKoleva et all (|2009l ) . On 
the other hand, their ages span a wide range. It is apparent 
that all dEs with blue central colours (dE(bc), blue symbols) 
are located towards fairly young ages, falling above the 3 
Gyr line in the figure. There is almost no overlap between 
the dE(bc) and the dE(N) (red symbols), except for VCC 
1353. Furthermore, most of the discy dEs (dE(di), green in 
colour) lie in between dE(N) and dE(bc). Thus, a clear di- 
vision of the different subclasses in terms of age is present, 
with low, intermediate and high ages for dE(bc), dE(di) and 
dE(N), respectively. Despite of having a much larger spread, 
the dE(N)s are clustered in the bottom left corner of the grid 
- or even below the grid in some cases - centered on the low 
metallicity region. Interestingly, only three dE(di)s fall left 
of the metallicity line of [Z/H] = —0.36 dex. 

Fig. HH presents the Mgb vs <Fe> diagram, where the 
data are superimposed on the TMB03 model with various 
a— element abundance ratios: [Q/Fe]= —0.3, 0.0, and 4-0.3 
dex. The results from previous studies of the stell ar popu- 
lations of low-mass system s ((Gorgas et al. 1997 : Gcha et al.l 
l2003l : iThomas et ai]|2003l ). that the a— element abundance 
of dEs is consistent with solar, seem to be in agreement with 
our study. However, the scatter is fairly large for the dE(N)s, 
and on average they might be slightly more q— enhanced 
than the dE(di)s. 

6.2 Age, metallicity and abundance ratios 

Using the combined sample of our galaxies and the re- 
measurement of the M08 galaxies' spectra with the method 
of the x^— minimization, we find a mean age of 3.0 ±0.8 
Gyr and a mean metallicity of [Z/H] = —0.31 ±0.10 dex for 
the dE(di)s. In contrast, the mean values for the dE(N)s are 
7.5 ±1.9 Gyr and [Z/H] = -0.54 ±0.14 dex. The dE(bc)s 
have comparatively young ages (i.e less than 3 Gyr). Two 
of them show a rather low metallicity (—0.90 ±0.25 dex for 
VCC0021 and -0.75 ±0.15 dex for VCC1912), while the 
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Figure 12. Central Mg versus <Fc> indices after correcting back 
to Lick resolution, superimposed on the TMB03 model with dif- 
ferent ages, metallicitics, and abundance ratios. The symbols are 
the same as in Fig. 1111 



third one also belongs to the discy subclass (VCC 0308), 
and indeed has a significantly larger metallicity (i.e —0.09 
±0.10 dex), consistent with the other dE(di)s. 

As a statistical comparison, we have used the 
Kolmogorov-Smirnov test for the goodness of fit with the 
null hypothesis that two observed distributions are from 
the same continuous distribution. The estimated probability 
Pxs that the difference is at least as large as observed if they 
had been drawn from the same population, is Pks = 0.01, 
0.03 and 0.15 for the age, metallicity, and [a/Fe]-abundance 
distributions, respectively, of dE(N)s and dE(di)s. This ap- 
parently confirms that the age and metallicity differ signifi- 
cantly for dEs with and without discs. However, this inter- 
pretation is challenged when the stellar population charac- 
teristics are considered with respect to galaxy luminosity; 
details are given at the end of this section. 

The galaxies in our sample follow a general trend link- 
ing age and metallicity: as displayed in Figure 1131 we see 
that the derived ages are decreasing with incre asing metal- 
licity of the dEs, similar to the observation of iRakos et al.l 
I2OOII for Fornax dEs. Interestingly, we see that an age- 
metallicity anti-correlation is pre sent for the dE(N)s, e ven 
tighter than previously noticed bv lPoggianti et al] (|200ll l for 
Coma cluster dEs. Such an age-metallicity anti-correlation 
is also found for giant early-type galaxies, but it is under 
debate because of the correlation of the metallicity and age 
errors (jTrager et al.l I2OO0I : iKuntschner et all 120011 ) . Never- 
theless, together with the dE(di)s, it can be noticed that the 
brighter dEs might even follow an opposite trend of increas- 
ing metallicity with increasing age. It needs to be mentioned, 
however, that the scatter of values is fairly large. 

In addition, we provide a plot with the error ellipses 
(Fig |14|l . to better illustrate the behaviour of our mea- 
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Table 3. Derived stellar population parameters using the model 
TMB03 



Galaxy 


Age, 


Gyr 


[Z/H] 


dex 


[a/Fe] 


dex 


VCC0308 


01.7 


+0.2 
-0.1 


-0.09 


± 


0.10 


0.01 




0.07 


VCC0216 


01.7 


+0.6 
-0.1 


-0.53 


± 


0.17 


0.19 


± 


0.13 


VCC0389 


04.0 


+0.8 
-0.8 


-0.21 


± 


0.07 


-0.01 


± 


0.08 


VCC0490 


02.7 


+0.9 
-0.6 


-0.16 


± 


0.12 


-0.04 


± 


0.11 


VCC0856 


03.5 


+0.9 
-0.4 


-0.21 


± 


0.07 


-0.01 


± 


0.08 


VCC0990 


03.6 


+0.3 
-0.6 


-0.19 


± 


0.07 


-0.06 


± 


0.08 


VCC1304 


04.1 


+ 1.1 
-1.1 


-0.65 


± 


0.17 


-0.04 




0.18 


VCC2019 


03.0 


+0.8 
-0.5 


-0.19 


± 


0.10 


-0.12 




0.11 



VCC0545 


11 


9 


+0.6 
-1.5 


-0 


97 


± 





05 


-0 


06 


± 





22 


VCC0929 


04 


3 


+0.9 
-0.4 





01 


± 





05 





00 


± 





05 


VCC1167 


10 


9 


+2.2 
-1.4 


-1 


00 


± 





07 





07 


± 





16 


VCC1185 


09 


9 


+3.1 
-1.3 


-0 


80 


± 





12 


-0 


25 


± 





16 


VCC1254 


06 


6 


+ 1.3 
-1.3 


-0 


41 


± 





12 





07 







08 


VCG1261 


02 


5 


+0.4 
-0.3 


-0 


07 


± 





07 





01 


± 





05 


VCC1308 


03 


8 


+ 1.2 
-0.6 


-0 


21 


± 





10 





07 







09 


VCG1333 


13 


7 


+ 1.3 
-5.4 


-1 


12 


± 





10 





12 


± 





23 


VCC1348 


15 





+0.0 
-0.7 


-0 


58 


± 





05 





38 







09 


VCC1353 


02 


4 


+0.9 
-0.7 


-0 


63 


± 





20 





19 


± 





22 


VCC1355 


03 


8 


+ 1.7 
-0.9 


-0 


43 


± 





22 


-0 


16 







16 


VCG1389 


10 


9 


+2.2 
-1.4 


-1 


02 


± 





07 





08 







18 


VCC1407 


14 


3 


+0.7 
-6.1 


-0 


70 


± 





15 





13 







13 


VCC1661 


08 


6 


+6.4 
-2.1 


-0 


73 


± 





20 


-0 


23 


± 





16 


VCC1826 


07 


5 


+ 1.9 
-1.8 


-0 


53 


± 





15 


-0 


15 


± 





14 


VCG1861 


03 


3 


+ 1.0 
-0.6 


-0 


09 


± 





10 


-0 


03 







08 


VCC1945 


03 


6 


+ 1.1 

-1 


-0 


41 


± 





17 





07 







12 



Galaxies from MichielsenetaL (2008 



VCG0021 


01.0 


+0.8 
-0.0 


-0.90 


± 


0.25 


0.00 




0.30 


VCC1912 


02.2 


+0.6 
-0.7 


-0.75 


± 


0.15 


-0.11 


± 


0.19 


VCC0397 


01.7 


+0.6 
-0.4 


-0.12 


± 


0.25 


-0.04 


± 


0.16 


VCC0523 


04.1 


+2.7 
-1.4 


-0.29 


± 


0.25 


0.15 




0.18 


VCC1183 


03.6 


+0.9 
-0.9 


-0.26 


± 


0.12 


-0.06 


± 


0.11 


VCC1695 


02.1 


+ 1.5 
-0.5 


-0.38 


± 


0.30 


0.04 


± 


0.26 


VCG1910 


05.7 


+2.5 
-1.4 


-0.09 


± 


0.10 


0.05 


± 


0.08 


VCC1947 


04.5 


+ 1.2 
-0.8 


0.13 


± 


0.10 


0.09 


± 


0.07 


VCC1087 


03.6 


+ 1.6 
-1.2 


-0.24 


± 


0.17 


-0.12 




0.16 


VCC1122 


05.2 


+ 1.7 
-1.9 


-0.61 


± 


0.25 


-0.29 


± 


0.20 


VCC1431 


13.1 


+ 1.9 
-2.7 


-0.53 


± 


0.12 


0.28 




0.12 


VCC1549 


04.3 


+ 1.4 
-1.5 


-0.12 


± 


0.15 


0.19 


± 


0.11 



surements with respect to the ag e and metallicity anti- 
correlation. iKuntschner et al]|200ll suggest that, when the 
measurement errors in the line strength are not neghgible, 
the tilt of the model grid leads to the correlated errors in 
age and metallicity, as visualized by the error contours in 
the age metallicity plane in Fig |14l 

The luminosity-weighted metallicities derived from the 
well measured nine indices are shown as a function of galaxy 
absolute magnitude in Figure lTSl (bottom left). Both dE(di)s 
and d E(N)s follow the well kn own metallicity— luminosity re- 
lation (|Poggianti et al.ll200ll ). In contrast, we do not see any 
strong correlation between measured metallicity and local 
projected galaxy density (bottom right). In the top panel of 
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Figure 13. Age versus metallicity [Z/H]. Here colour and sym- 
bols are the same as in Fig llOl 



the Figure, one can see that the dEs in general are consistent 
with solar [a/Fe]-abundance, and we do not see any relation 
with luminosity or local density. 

Unlike the correlation of metallicity and luminosity, the 
measured ages show an aniti-correlation with the luminos- 
ity, and correlation with local projected density, although 
the latter trend is not as strong as the former. It appears 
that this might be the cause of the apparent age difference 
between dEs with and without discs, since in addition to 
the relation with magnitude, the dE(N)s are situated in the 
denser cluster regions as compared to the dE(di)s. Thus, 
taking the left and right panels in Fig. [15] together, no sig- 
nificant difference is seen between dE(di)s and dE(N)s in 
either age or metallicity, as long as they are compared at the 
same magnitude and local density. 

On the other hand, the bimodality in the age distribu- 
tion is quite prominent. Interestingly, the majority of those 
dEs belonging to the group with lower ages lie consistently 
within the one-sigma area of the dE(di)s (gray-shaded region 
in Fig. 15). But we do not see any such strong bimodality 
in the metallicity distribution. 



7 DISCUSSION AND CONCLUSION 

We have undertaken a spectroscopic investigation of Virgo 
cluster dEs, deriving ages, metallicities and a— element 
abundance ratios for 38 dEs of different types, i.e., dEs with 
and without discs, and examined the variation of their stellar 
populations according to their morphology. We have estab- 
lished that the mean ages and metallicities of different types 
of dEs are different. 

This result mostly represents the central part of the 
galaxy, as we used a central fixed aperture to extract the 
one dimensional spectra of the galaxies, in order to optimize 
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Figure 15. The derived ages (top), motallicities (middle) and [a/Fe]-abundance (bottom), plotted against r band absolute magnitude 
(left), and against local projected density (right, number per sq.deg.). Colours and symbols arc the same as in Fig llOl The gray shaded 
region represents the mean iter age and metallicity of dEs(di) 



the SNR for our Lick index analysis. We expect tiiat having a 
wider aperture effectively minimizes the influence of central 
nucleus. In a worst-case scenario, for those dEs that have a 
relatively low surface brightness, i.e., a rather faint stellar 
envelope compared to the nucleus, it is not unquestionable 
whether our extracted spectra do really represent the galax- 
ies as a whole, or the nuclei only. Interestingly, though, we 
found that most of the faint nucleated dEs have a fairly 
old stellar population, and that typi cally, the nuclei tend t o 
be younger than their host galaxies jPaudel fc Liskei 2009h. 



This appears consistent with the stud i es of IChilingarianI 
(|2009l ') for Virgo dEs and iKoleva et all (|200g ) for Fornax 
cluster dEs: both concluded that the central part of the 
galaxies are relatively younger and more metal-rich than the 
outer part. In that respect, we can argue that our conclu- 
sions are not biased by having a strong effect of the central 
nucleus in the spectral extraction process. If anything, our 
result of old galaxy ages at fainter magnitudes could only be- 
come somewhat more pronounced by removing the nucleus 
contribution. 

It was known before that the different dE subclasses, 
in particular those with and without discs, exhibit differ- 
ences in their colours, magnitude, and the local environment 



in which they reside (|Lisker et al.ll2007h . The dE(di)s were 
found to ha ve, on averag e, slightly bluer colours than the 
dE(N)s, and iLisker et al.l interpreted this colour difference 
with either older age or higher metallicity of the dE(N)s. By 
combining our sample with that of M08, and therefore reach- 
ing statistically significant subsample sizes, we see that the 
colour offset is mainly an effect of age rather than metallic- 
ity. However, it is not straightforward to tell which govern- 
ing factor is most responsible for this apparent correlation 
of galaxy substructure and stellar population characteris- 
tics. As noted above, we find correlations of age with both 
magnitude and local density, and we know that the dE(di)s 
are relatively bright and typically located in less dense clus- 
ter regions than the dE(N)s; the combination thus shows 
that no significant difference between the stellar popula- 
tion of dE( di)s and dE(N)s ca n be claimed. This is inter- 
esting, since [Lisker et ahl (|2008l ) pointed out that the fairly 
wide (estimated) distribution of intrinsic axial ratios of the 
dE(N)s would be consistent with assuming that the dE(di)s 
- most of which are nucleated - are simply the flat tail of 
this distribution. The location in regions of comparatively 
lower density could then be interpreted such that discs in 
dwarfs can only survive for a reasonably long time if the 
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Figure 14. Error contours in different projection planes of age, 
metallicity and [o/FeJ-abundance space. 

am ount of dynamical heatin g due to tidal forces is low (also 
see iMastropietro et al-llioOSl ). and therefore, dEs with discs 
are preferentially found outside of the cluster center. More- 
over, L06 noted that there could still be a significant fraction 
of objects that have discs, but that were simply not found 
by their analysis of SDSS images, due to low SNR or the 
lack of disc features like spiral arms or bars. We see in our 
data that a handful of dE(N)s are as young and metal rich as 
the dE(di)s, and interestingly, they also lie in a low density 
region. So maybe, just as a possibility, these could be those 
dE(N)s that actually host discs, but were not identified as 
such — which would then also explain the bimodality when 
compared to the fainter ones. 

Furthermore, an apparent inconsistency is observed in 
the relation between age and luminosity: previous studies 
provided evidence for a continuously later formation epoch 
and/or a longer duration of star formation w hen going from 
the massive ellipticals to the dwarfs (see, e.g. , iGavazzi et al.l 
I2OO2I : iBoselh et al.|[2005l : iLisker fc Hanll2008h . which appar- 
ently is at variance with our findings. On the other hand, 
this trend is unlikely to continue to even fainter magnitudes, 
since most of the Local Group dwarf spheroidal (dSph) 
galaxies are know to be dominated by old populations (e.g. 
iGrebelll 19991 ). It seems that, the fainter dEs are significantly 
older than the brighter ones, no matter at which density 
they are located. To statistically confirm this result, we di- 



vide the total sample in two groups: those being fainter and 
brighter than the mean magnitude M,- = —17.22 mag. We 
obtain a mean age of 3.9 Gyr and 7.4 Gyr for the brighter 
and fainter group, respectively, which is very similar to the 
mean ages of the discy and nucleated subsamples (3.0 and 
7.5 Gyr). Furthermore, the K-S test supports our finding 
that a dichotomy in the distribution of ages can be seen not 
only with respect to different morphology but also different 
luminosity, as we get a K-S test probability of 0.01 for the 
null hypothesis that the age distributions of the fainter and 
brighter group are not different. 

The ages of 10 Gyr and above indeed appear more com- 
parable to the Local Group dSphs. In that sense, we find a 
discontinuity in the trend of stellar population characteris- 
tics with luminosity, and we can speculate that the fainter 
dEs might be a different species than the brighter ones or 
simply were not able to sustain star formation for as long a 
period as the more luminous and presumably more massive 
dEs. 

In comparison to previous st ellar population studies o f 
dEs in different clusters, unlike IChilingarian et al] (|2008l ) 
who found a correlation of [a/Fe] with the projected distance 
in the Abell496 cluster - although o nly in the very central 
part - and the study of lSmith et al.l (|2009l ) for Coma, we do 
not find any correlation between a-element abundance ratio 
with environmental density. This could also be due to differ- 
ences between the clusters, in the sense that Coma is a more 
virialized cluster of higher richness class than Virgo. Never- 
theless, our range of age and met allicity roug h ly agr ees with 
their ranges. Another study by ICeha et al.l (|2003l ) , focus- 
ing on the stellar populations of rotating and non-rotating 
Virgo dEs, resulted in no difference between them. On the 
other hand, a recent, more extensive study bv lToloba et al.l 
liooi) finds that rotating Virgo dEs are typically located 
further from the cluster center and have younger stellar 
population ages. We have only six dEs in common with 
these studies (counting only those for which tabulated val- 
ues were published). Three of them (VCC 0308, VCC 0856 
& VCC 1947) have disc substructure and are rotating, while 
for three (VCC 1254, VCC 1261 & VCC 1308) no disc fea- 
tures were i dentified and they are non-rotating according to 
ICeha et al.l . Note, though, that VCC 1261 is clearly rotat- 
ing fro m the analysis of dE globular clusters of lBeaslev et al.l 
( 2OO9I) . [ChilingarianI (|2009l ) shows that this galaxy contains 
a rotating kinematically decoupled core, indicating the diffi- 
culties in drawing robust conclusions on the global rotation 
of a dE from the inner stellar kinematics. We also find that 
the overall stellar populations are quite mixe d within these 
six ga laxies, as has already been noted by Ivan Zee et al.l 
(|2004h in their study of (non-)rotating dEs. 

There have been numerous discussions on the origin of 
the passive dwarf galaxies. The key issue is whether cluster 
dwarfs are primordial, having formed from the h ighest den- 
sity peaks in the proto-cluster (|Tu11v et al.ll200 and avoid- 
ing subsequent merging, or originating from transformation 
of late-type disc galaxie s through the interaction with the 
cluster environment (e.g. lBoselli fc Gavazzi,2006 V However, 
to explain the formation of dEs, using such single scenarios 
may be too simplistic, and it is likely that both scenarios 
operate at some level. Furthermore, the question is how to 
find unique features for such a process that dominates in 
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a given mass regime, a given environment, and for a given 
morphological subclass of dwarfs. 

Our results promote the idea of different scenarios for 
different morphology and/or different luminosity of dEs, al- 
though the physical mechanisms responsible for the quench- 
ing of the star formation activity in the galaxies are not 
easy to understand from the stellar population properties. 
The observed systematic difference in age & metallicity be- 
tween different types of dEs support the finding of L07. The 
typical age range of the brighter dEs of less than 5 Gyr in- 
dicates that the star formation activity ceased at z < 0.5. 
The younger age, higher metallicity and consistency with 
solar Q— element abundance would suggest that a discy dE 
can indeed form through the structural transformation of a 
late-type spiral into a spheroidal system, triggered by the 
popular scenario of strong tidal interactions with massive 
cluster galaxies. A stellar disc component (which has rather 
high metallicity) may survive and form a bar and spiral fea- 
tures that can be retained fo r some time depending o n the 
tid al heating of the gala xy (|Mastropietro et al.l l2005l : also 
see iLisker fc Fu"ch3l2009l ). This "gal axy harassment" pro- 
cess could thus produce disc-shaped dEs. 

On the other hand, the fainter dE(N)s, with their older 
ages, might be a different class of objects with a different 
formation scenario. The poorer metal content supports the 
ide a that they might be prim ordial objects, as suggested 
bv lRakos fc Schomberj (|2004f ): they might have either suf- 
fered early infall into the cluster potential, or formed to- 
gether with cluster itself. The generally low surface bright- 
ness of dEs suggests that whichever mechanism is responsi- 
ble for the halting the star formation activity, it must have 
been rather efficient. The common idea is that internal feed- 
back might be responsible for their cessation of star forma- 
tion activity at such early epochs. We can speculate that 
the slightly higher mean [a/Fe] as compared to the dE(di)s 
could indicate a more efficient star formation activity than 
their discy counterparts before the quenching. High-efiiency 
star formation has also been suggested as an explanation for 
the comparatively high metallicity of the old populations of 
Local group dSph as compared to dirr populations of the 
same age, i.e., as an explanatio n for the offset alo ng the 
metallicity— luminosity relation (jCrebel et al.l [20031 ). How- 
ever, further dynamical studies of these fainter galaxies 
are needed to prove whether the removal of gas by inter- 
nal feedback processes is efficient or not, since it has been 
pointed out that it is difficult to really eject the gas from 
a dwarf galaxy by supernovae explosio ns unless the dwarfs 
have masses less than lO''' solar masses (|Mac Low fc Ferraral 
I1999D . 
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APPENDIX A: LICK OFFSETS 

For this project and instrument configuration there were 
no stars obser ved in common with the original Lick/IDS 
stellar library (|Worthev et al.lll994bl ) which is the standard 
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Figure Al. Comparison of Lick/IDS index measurements with the MILES stellar library (Sanchez06) for stars in common (plus signs). 
The red dashed line in each panel shows the mean offset derived by a biweight estimator (see also Table EU l . In the title of each panel 
we list the index name, the index offset and associated error, a robust estimate of the la scatter and the number of stars used in the 
comparison. A subset of the indices show evidence for line-strength dependent offset. For indices where this is a significant effect (2(t) 
we also show a robust fit to the relation (green dash-dotted line). See text for details. 
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way to calibrate the inde x measurements to the Lick system 
by deriving small offsets (|Worthev fc Ottavianill 19971 '). How- 
ever, in t his paper wc made use of the MILES stellar library 
ISanchez-Blazgucz ct al. 200^ ) to determine the offsets for 
indices measured on flux calibrated spectra to the original 
Lick/IDS system assuming that the MILES library and our 
data are well flux-calibrated. 

The spectra were flrst broadened to the Lick/IDS res- 
olution with a wavelength-dependent Gaussian assuming a 
spectral resolution of the MILES library of 2.3 A (FWHM). 
Then, we compared the index measurements in common 
with the Lick stellar library for up to 239 stars depending 
on the availability of the measurements. The offsets and as- 
sociated errors (see Figure lAll and Table IA1[) were derived 
with an outlier resistent biweight estimator (IDL Astron- 
omy library biweightjnean.pro). The offsets are generally 
small (<0.2 A) but individual indices show larger offsets 
(e.g., Fe4383, Fe5015). Additionally, several indices show ev- 
idence for line-strength dependent offsets (e.g., CN2, G4300, 
Mg2; see also Vazdekis 1999). We used an outlier-resistant 
two-variable linear regression method (IDL Astronomy li- 
brary robust_linef it .pro) to investigate the significance 
of such trends. For indices where we find a slope different 
from zero with more than 2a signiflcance we show the flt 
in Figure lAll with a green dash-dotted line and give the flt 
parameters in Table [AT] 

We note th at previous attempts to es tablish Lick 
offsets with the JonesI 1997 stellar library (jNorris et al.l 
l2006llWorthev fc Ottavianil 1 19971 ) have resulted in signifl- 



Table CI. Derived stellar population parameters using the 
AV/MILES model at 11 A. 



cantly different offsets for some indices, perhaps most no- 
tably the higher order Balmer lines. We ascribe the differ- 
ences in offsets t o an imperfect flux calibration of the Jones 
library. Recentlv. lSanchez-Blazguez et alll2009l also used the 
MILES library to derive Lick offsets and flnd similar results. 
However, signiflcant differences remain for individual indices 
(e.g., Mg6) where we favour our outlier robust analysis. 

For the present study, we apply the offsets listed in Ta- 
ble [AT] (column 2). Only one of the indices (Fe5015), used 
in this study to derive the stellar population parameter esti- 
mates, shows evidence for a line-strength dependent offset. 
Since the trend is very weak in the range of observed line- 
strength for this index we apply only a normal offset. 
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APPENDIX B: LICK/IDS INDICES 

In Table [BT] we list the measured indices and their errors, 
which we used for the SSP parameter estimation of the 
galaxies of our sample. The indices are transformed to the 
Lick/IDS system using the offsets described in Appendix [A] 



APPENDIX C: COMPARISON OF THE 
ESTIMATED SSP PARAMETERS 

In this section we present several comparisons of the SSP 
parameters derived from different SSP models, using differ- 
ent methods and also different data sets. We also provide 
the derived age and metallicity from the AV/MILES model 
for our data set in Table [CT] 
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Table Al. Offsets to the Lick/IDS stellar library system derived from the MILES library. 
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Note. Column (1) gives the index name, while column (2) gives the mean offset (Lick - MILES) to the Lick/IDS 
system. Column (3) gives a robust estimate of the Itr scatter and column (4) shows th number of stars used in the 
comparison for each index. Columns (5) and (6) give the coefScents of a linear fit to (Lick - MILES) = A + Bx Lick 
where a significant slope was found. 

Table Bl. Measured central line-strength indices corrected to the Lick/IDS system, measured in the central 4 arcsec with one sigma error 
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Figure C2. Comparison between the Lick indices measured from our and the M08 sample for the galaxies common to both samples, at 
the resolution 11 A 
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Figure C3. Lick/IDS index— index diagrams of the age-sensitive Balmcr indices H/3, Hfp and H5p versus the metallicity sensitive 
indices Mgb,<Fe> and [MgFe]', with overplotted model grids of TMB03 
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Figure CI. Comparison between the derived age (left) and 
metallicity (right) with TMB03 for our sa mple and that o f M08 . 
We use different approa ches, i.e. following iProctor et al. I ||2004| ') 
and lCardiel et aP l|2003h by us and by M08, respectively. 
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Figure C4. Comparison of derived stellar population parameters 
from the AV/MILES model at different resolutions for the M08 
sample. 



